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Abstract 

A constraint between thermal fireball parameters arises from the requirement that 
the balance of strangeness in a fireball is (nearly) zero. We study the impact of this 
constraint on (multi-)strange (anti-)baryon multiplicities and compare the hadron gas 
and quark-gluon plasma predictions. We explore the relation between the entropy 
content and particle multiplicities and show that the data are compatible with the 
quark-gluon plasma hypothesis, but appear to be inconsistent with the picture of an 
equilibrated hadron gas fireball. We consider the implications of the results on the 
dynamics of evolution and decay of the particle source. 
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1 Introduction 



The reported enhancement of strange particle production in heavy ion collisions relative 
to proton-proton and proton-nucleus collisions |1[ tells us that we are dealing with a state 
of matter which is i) very dense or ii) relatively long-lived, or in which in) strangeness 
production cross sections are enhanced, or which iv) possesses some combination of these 
three factors. In order to be able to be more specific about the nature of the dense matter, 
the relative production abundances of strange and multi-strange baryons and anti-baryons 
were studied. They turned out to be particularly sensitive probes of the thermal conditions 
of their source [§, g|. Studies of the properties of a hadron resonance gas fireball including 
strange particles have shown ||, |, |6], |7|] that the constraint of vanishing net strangeness of 
the emitting source leads to similar thermal and chemical conditions for both a hadronic gas 
(HG) and a deconfined quark-gluon plasma (QGP) fireball, if the temperature is taken to be 
near T = 210 ± 10 MeV. This value was extracted || from the m_|_-spectra from 200 GeVA 
S-W and S-Pb collisions under the assumption that their slope is dominated by the thermal 
motion of the particles at freeze-out and does not contain a collective flow component ||. 
We will discuss that the observed smallness of the strange quark chemical potential appears 
to be inconsistent with the presence of a strong flow component in the spectra. 

It has already been pointed out |J that near T = 210 MeV further information, such as 
the particle multiplicity per participant nucleon, may help to distinguish QGP from the HG. 
Following this route, we will show here in detail [JTUJ] that the data of emulsion experiment 



EMU05 [FT] indeed appear to be incompatible on these grounds with the hadronic gas 
fireball model. Further support for a new physical state such as the quark-gluon plasma as 
the source of these anti-baryons may be obtained by studying the response of the measured 
parameters to changes in beam energy or size of the colliding nuclei. Our work has been 
in part motivated by the imminent extension of the experimental program at CERN and 
BNL to heavy nuclear projectiles (Pb and Au, respectively), and by the prospect of higher 
(and lower) beam energies becoming available, in particular towards the end of the decade at 
Brookhaven with RHIC and possibly at CERN with LHC. Therefore, in addition to T ~ 200 
MeV appropriate for the CERN-SPS energy range, we also include in our discussion the 
parameter range (T ~ 150 MeV) suitable for the current BNL heavy-ion experiments at 
10-15 GeVA and for potential future low beam energy runs at CERN, as well as T ~ 300 



MeV, which we judge appropriate for the future RHIC and LHC facilities JL2 



Our expectation is that, as the conditions of the experiments are changed, we will be 
able to study the response of the statistical parameters such as the strange quark chemical 
potential /i s or the strangeness saturation factor 7 S . These parameters show characteristic 
differences in a hadron gas and a QGP: 7 S , for example, is expected to rapidly approach unity 



for increasing four-volume of the QGP fireball, while remaining small in HG fireballs |jl3 



The strange quark chemical potential /i s , on the other hand, has the characteristic property 
that it is exactly zero in a strangeness-zero QGP fireball, independent of its baryon density. 
In contrast to this, /j, s is generally different from zero in conventional hadronic gas reactions, 
except for the special case of baryon-free matter with /ib = 0; however, if T is not too large 
(T < 220 MeV), for each value of the temperature a certain value of baryochemical potential 
//b 7^ can be found for which also /if G = 0. This value generally depends on the detailed 
form of the equation of state (EoS), and for the case of a conventional (Hagedorn type) HG 
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it has been discussed in Refs. [0, |H| . While the large values of 7 S which are characteristic 
of the fast chemical kinetics in a QGP are largely preserved during hadronization, the value 
of n s is usually subject to modifications, in particular if the phase transformation occurs 
slowly. Experimentally, the values of /i s and 7 S are reflected in the observed strange hadron 
ratios (where, as we will discuss, in the case of QGP certain model assumptions for the 
hadronization process are necessary in order to quantify this relationship), which allows for 
a rather straightforward experimental test of the fireball picture, of the concept of (absolute 
or relative) chemical equilibrium, and of the nature (HG or QGP) of the collision fireball. 

In order to distinguish between the two phases by studying particle production as the 
conditions of the experiments are changed we need to understand also in the HG phase 
how the strange chemical potential /i s relates to the baryochemical potential /^b in the 
collision fireball. Such a relationship is required by the fact that the total strangeness in the 
initial state is zero, and that strangeness is conserved by strong interactions. Thus, up to 
(presumably small) pre-equilibrium emission effects which can introduce a small strangeness 
asymmetry s — s = es, the number of strange and anti-strange quarks in the fireball are 
equal (e ~ 0). We work out in detail how the strange chemical potential is determined 
at fixed temperature by the baryochemical potential and a given (small) fraction e of net 
strangeness in the fireball. We find that at a fireball temperature around T = 210 MeV, 
there is a peculiar behavior of fi s which causes a particular insensitivity of the strange anti- 
baryon ratios alone to the structure of the source. In order to be able to reach a conclusion 
about the internal structure of the source, we then study the entropy content of the particle 
source. We also show that the kaon to hyperon ratio is a poor signature for the nature of 
the source, but could be a good measure of the baryochemical potential |~5] . 



The HG part of our work is similar in spirit, but different in detail from the parallel effort 
of Cleymans and Satz [^]0. We study particle ratios at fixed transverse mass m± > 1.7 GeV. 
We allow for the strangeness phase space to be only partially saturated and show how the 
degree of saturation 7 S can be experimentally measured in order to later serve as a test for the 
kinetic theory of strangeness production. We allow the strangeness to be slightly unbalanced 
(up to 10%) and discuss the corresponding variations in the observables. Furthermore we 
compare the entropy content for the two different scenarios (HG and QGP) and confront the 
findings with the observed particle multiplicity. In particular we consider the ratio of net 
charge to total charged multiplicity, Dq = (N + — N~)/(N + + N~), which can be measured 
(without identifying particles) with any tracking device within a magnetic field; we show 
that this ratio is directly related to the specific entropy S/E of the fireball, which should be 
considerably larger for a QGP than for a HG 0. Since, due to the different stopping, the 
experimental results from S-S and S-W collisions at 200 GeV A are not in the same class, we 
avoid combining these results because this would weaken the conclusions from our analysis. 
We do not combine in our analysis the global kaon data with anti-hyperon data, since low p± 
kaons, unlike strange anti-baryons, can also arise from peripheral, spectator related processes 
and other rescattering processes and are therefore not necessarily messengers from the same 
stages of the collision. Instead we consider the A/p and A/p as well as and Q/E~ 



1 When comparing the present work with Ref. Q], it should be noted that these authors use as the 
strangeness potential the chemical potential /ig of the (anti-)hyperon number, while we prefer to work 
directly with the strange quark chemical potential /i s . The relationship between their /ig and our /i s is given 
by fi s = /i B /3 - /^s- 
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ratios which, as we will explain, provide independent tests of the degree of strangeness 
saturation 7 S . 

At this point it is worth recording that if a QGP state was formed (as may be indicated by 
the relatively large value of 7 S , see below), the strangeness pair abundance per participating 
baryon in this state would be about a factor 2.5 larger compared to hadronic gas interactions. 
On the other hand, there would be a similar enhancement in the produced entropy (which 
will manifest itself in the total multiplicity), resulting in a difficulty to identify the absolute 
enhancement of both quantities: the simplest observable, the K/tt ratio, is therefore not a 
very good discriminator in this context ||16|| . As we will show, measurement of the specific 



entropy S/B, either globally via the total particle multiplicity per participating nucleon, or 
(better, see below) in the same kinematic range where the produced strange particles are 
observed via the ratio Dq = (N + — N~)/(N + + N~), can remove this ambiguity. 

One of the important objectives we pursue is to present a large number of predictions for 
(strange) particle production which are capable to test the internal consistency of the model 
developed here and can shed new light on the problem of distinguishing the HG and QGP 
fireballs. These predictions are quantitative for the S-W collision system at 200 GeV A, but 
are of a more qualitative nature for the yet unexplored Pb-Pb collisions where we cannot 
predict precisely the governing statistical parameters (temperature, chemical potentials). We 
also discuss the systematic behavior of the different statistical model observables in changing 
environments. 

Perhaps the most notable conclusion we reach concerns the model of hadronization of 
the fireball: we conclude that the only consistent model of the reaction involves explosive 
disintegration of a QGP as the source of the observed strange anti-baryons. This allows the 
chemical potentials of the quarks in the QGP phase to get transferred without modification to 
the observed high-m^ strange baryons and anti-baryons, via quark coalescense, and explains 
naturally the observed vanishing of the strange chemical potential. Since only high m± 
domain of strange anti-baryon spectra has been experimentally explored, the simple picture 
of a complete explosion is not absolutely necessary. It is conceivable that at low m±, in 
particular for low-p^ pions, another slower hadronization mechanism also plays a role. We 
do see the need for some mechanism that produces additional pions, in excess over those 
obtained from the explosive quark coalescence picture, in order to save the entropy balance 
and to account for the relatively large observed multiplicity densities. These excess pions 
(plus perhaps some other light mesons) could easily arise from the gluons in the QGP which 
carry a large fraction of its entropy. A quantitative formulation of such a hadronization 
model, however, is not yet available and would require additional data on the detailed origin 
of the excess multiplicity and its distribution in momentum space. 

Naturally, these conclusions rest on the assumed validity of the thermal fireball model. 
While at present the number of parameters just equals the number of independent observ- 
ables, the model has considerable predictive power and, with the help of the predictions 
presented in this paper, its validity will be (hopefully soon) tested. Other pictures for the 
source of the observed hadrons, such as an equilibrated HG or a QGP fireball which evolves 
slowly and in full equilibrium into a HG fireball, can not account consistently for certain 
important experimental features of strange anti-baryon production. 

In section || we present our thermal model with special attention given to the parameters 
describing the chemical equilibrium. We take into account the u-d asymmetry and study the 
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condition of stangeness balance. In section [5] we determine the thermodynamical parameters 
using the WA85 data on strange baryon and anti-baryon production (with and without 
taking into account resonance decays). In section |] the implication of these parameters on 
the entropy content of the QGP and HG fireballs and on the multiplicity are discussed. 
We then confront our model with EMU05 and NA35 data. In section [| we discuss the 
implications of our results regarding the nature of the initial fireball and the dynamic of its 
evolution. We give some concluding remarks concerning the distinction we draw between 
QGP and HG. 



2 Thermal models 

The use of thermal models to interpret data on particle abundances and spectra from nuclear 
collisions is motivated by the hope and expectation that in collisions between sufficiently 
large nuclei at sufficiently high energies a state of excited nuclear matter close to local 
thermodynamic equilibrium can be formed, allowing us to study the thermodynamics of 
QCD and the possible phase transition from a HG to a QGP at a critical energy density. 
Since we do not yet reliably know how big the collision system and how large the beam 
energy has to be for this to occur (if at all), thermal models should be considered as a 
phenomenological tool to test for such a behavior. Even if successful, the validity of such 
an approach must be checked later by a more detailed theoretical analysis of the kinetic 
evolution of the collision fireball. 

It has been noted repeatedly that, especially in collisions involving heavy nuclear tar- 



gets, the observed particle spectra, after correcting for resonance decay effects jl7| , resemble 
thermal distributions with a common temperature T (although a flow component cannot 
presently be excluded and is actually quite likely ||), and that so far the interpretation in 
terms of a generalized thermal model, which assumes equilibration only of the momentum 
distributions and the light quark abundances, but allows for deviations in particle abun- 
dances from absolute chemical equilibrium, in particular in the strange sector, appears to be 
consistent with the observed particle ratios ||. In this section we will discuss this model as 
far as it is relevant for our analysis. 



2.1 Fireball parameters 

Let us assume that in collisions between two nuclei a region with nearly thermal equilibrium 
conditions is formed near "central" rapidity, i. e. at rest in the center of momentum frame 
of the projectile and a target tube with the diameter of the projectile. 

This central fireball from which the observed particles emerge is described by its tem- 
perature T and by the chemical potentials /ij of the different conserved quark flavors u, d, s. 
Since the conserved quantum numbers of hadrons are simply the sum of the corresponding 
quantum numbers of their quark constituents, the hadronic chemical potentials are given 
as the sum of constituent quark chemical potentials. Since the strong and electromagnetic 
interactions do not mix quark flavors, u, d, and s quarks are separately conserved on the 
time scale of hadronic collisions, and they can only be produced (or annihilated) in pairs: 
X <-> q^i. This implies that in absolute chemical equilibrium the chemical potentials for par- 
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ticle and anti-particle flavors are opposite to each other. This observation can be restated in 
terms of the particle fugacities which are convenient for counting particles and are related 
to the chemical potentials by: 

k = e^ T . (1) 

Since the quantum numbers of a hadron are simply the sum of the corresponding quantum 
numbers of their constituent quarks, the fugacity of each HG species is simply the product 
of the quark fugacities, viz. A p = A^Ad, Ak+ = A U A§, etc. Consequently the relation between 
particle and anti-particle fugacity is in general: 

Aj = A" 1 . (2) 

If absolute equilibrium is not reached, this relationship can still be maintained by in- 
troducing an additional parameter characterizing the approach to equilibrium. For strange 
particles, specifically, absolute chemical equilibrium is not easy to obtain: due to the higher 
mass threshold, the production of ss pairs usually proceeds much more slowly than that of 
uu and dd pairs. This is true both in hadronic interactions and in the quark-gluon plasma 
[O]. In fact, the ss equilibration time scale is estimated to be of the same order of magni- 



tude as the total lifetime of the collision fireball in the QGP phase [18], and much longer 



than that in the hadronic phase fl9|] . On the other hand, once a certain number of strange- 
anti-strange pairs has been created in such a hadron gas, their redistribution among the 
various strange hadronic species is no longer hindered by any thresholds and occurs rela- 
tively fast. It thus makes sense to introduce the concept of relative chemical equilibrium, in 
which the strange phase space is not fully saturated, but whatever strangeness has already 
been created is distributed among the available strange hadron channels according to the 
law of maximum entropy. Parametrizing the approach to strange phase space saturation by 
a factor < 7 S < 1, the individual hadron fugacities are then given by 

lt \ = Ai , (3) 

where the \ still satisfy the chemical equilibrium condition (0), and where the power n(i, s) 
counts the total number of strange valence quarks plus anti-quarks in hadron species i. 

Altogether, within the thermal model the fireball is thus described by 5 parameters: its 
temperature T, its baryon chemical potential 

A*b = 3/i q , (4) 
the strangeness chemical potential fx g , the light quark asymmetry potential 5/j, (see subsec- 



tion p.2|) and the strangeness saturation factor 7 S . The value of the temperature can be 
extracted from the observed transverse momentum spectra whose exponential slope is gov- 
erned by a combination of thermal motion and collective expansion of the fireball ||. The 
separation of these two contributions, i.e. the isolation of T, is not straightforward, but has 
been studied extensively in 0. Given T, the number of independent variables is further 
reduced by studying the u-d asymmetry in the fireball and relating it to the asymmetry 
potential 5{i, and by the condition of strangeness conservation. This leaves us with two 
independent variables, say /xb and 7 S , which have to be determined from measured particle 
ratios. 
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The various parameters described above enter the thermodynamic description through 
the partition function from which all observables are derived. For the hadronic gas state 
we follow here the conventional hadronic gas model first developed by Hagedorn [20| in 
which the equation of state of any hadronic system is obtained from a partition function 
which sums over all hadronic resonances. It is believed that by including a complete set 
of intermediate scattering states in the form of (usually zero-width) hadronic resonances 
the dominant interactions within the hadronic system are automatically taken into account. 
However, an important but tacit approximation is made in this approach, concerning the 
parameters such as mass (and sometimes widths) of the resonances which are inserted with 
their free space values. This neglects the possibility that in particular at high temperature 
and particle or baryon density, there may occur significant shifts from the free space values- 
one often speaks in this context of the melting of hadronic masses as one approaches extreme 
conditions. While we do our computations entirely within the conservative framework of the 
conventional model, our general physics discussion will be shaped in such a way that our 
conclusions should not critically depend on the exact validity of the model used. 

The standard hadronic gas model is certainly not sufficiently precise for our purpose if we 
were to ask questions concerning the absolute values of energy, entropy, baryon density, etc. 
which depend either on the ever increasing mass spectrum of particles or on the proper vol- 
ume occupied by the particles |2T| . However, the condition of zero strangeness and quantities 



such as the entropy per baryon, which involve ratios of extensive variables, are independent 
of the absolute normalization of the volume and of the renormalization introduced by the 
diverging particle spectrum and hence can safely be considered in our approach. 

We now turn to the discussion of the constraints which limit the number of the free 
thermal parameters in the fireball. 

2.2 u-d asymmetry 

In our calculations we will distinguish between the u and d quarks by introducing separate 

fugacities, A u and Ad, instead of A q in the partition function. Since even in the heaviest 

nuclei there is only a small u-d asymmetry due to the neutron excess, we find it convenient 
to work with the variables: 

Mq= (Hd + Hn)/2 , (5) 

<5/i = /i d - fi u ■ (6) 

Here, /i q is the "light quark" chemical potential (/ib/3 in terms of the baryochemical po- 
tential) and 5fi describes the asymmetry in the number of up and down quarks due to the 
neutron excess in the heavy ion collision. The ratio of the net number of down and up quarks 
in the fireball is 

{u — u) 

In a central S-W collisions, where a tube with the transverse area of the S projectile is swept 
out from the W target and participates in the fireball Rf~ w ~ 1.08; in Pb-Pb collisions 
.RPb-Pb = 1.15. The value of 5fi is at each fixed T determined by the value of Rf, but 
depends on the assumed structure of the source, i.e. the EoS. We have investigated this 
relation for the two cases considered here: the conventional HG and the QGP. 
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In the HG phase we include in the par- 
tition function all non-strange mesons up to 
a mass of 1690 MeV, all nucleons up to 1675 
MeV and all A's up to 1900 MeV. [We note 
that higher resonances would matter only if 
their number were divergent (as is the case 
in the Bootstrap approach of Hagedorn [BU| ) 
and if the HG was sufficiently long lived to 
populate all high mass resonances]. In ad- 
dition strange hadrons are included as de- 
scribed in Eq. ( |10"D below, which we expand 
in calculations to allow to distinguish be- 
tween u and d quarks within the various 
strange hadron multiplets. 

In Fig. p] we show how in the conven- 
tional hadron gas the difference Sfi/fiq in- 
duced by the isospin asymmetry of the fire- 
ball depends on the temperature, for two 
selected values Rf = 1.08 and 1.15 (appro- 
priate for S-W and Pb-Pb collisions, re- 
spectively), and for three fixed values of 

A s = 0.95, 1.00, and 1.05. The curves shown correspond to non-trivial solutions (i.e. solutions 
with /i q ^ 0) of the strangeness balance equation with the given values for A s ; as we will 
explain below, these solutions cease to exist (since we require (s — s) ~ 0) above a certain 
critical temperature, which is given by T ~ 230 MeV for A s = 1 and smaller (larger) for 
smaller (larger) values of A s . We see that, while for T — > the ratio S/j,//j, q vanishes in the 
hadron gas, it is of the same order as the deviation from unity of Rf once T reaches values 
around 200 MeV which are of interest here. 

For the QGP the ratio Sfi//i q is independent of A s , due to the decoupling of the strange 
and non-strange chemical potentials in the partition function. For /x q < tcT one finds [[| 



200 
T [MeV] 



Figure 1: 5fJ,/u. q asymmetry for Rf = 1.08 
(S-W collisions) and Rf = 1.15 (Pb-Pb colli- 
sions). Solid lines: A s = 1.00; upper dotted lines: 
A s = 0.95; lower dotted lines: A s = 1.05. 



QGP 



1 



5fi 



(8) 



such that 5(j,//j, q is exactly equal to -Rf — 1. It is remarkable that towards the maximum 
temperature consistent with strangeness conservation the HG and QGP based results for Sfj, 
at a given Rf nearly agree. We thus find that, although small, the difference between the 
chemical potentials of u and d quarks is not always negligible. We note for later applications 
that in the region of interest to us here (T ~ 210 MeV) we have irrespective of the nature 
of the fireball: 

SfJ. yUq 



T T 



(Ri ~ 1) • 



(9) 



2.3 The strangeness partition function 

In the Boltzmann approximation, it is easy to write down the partition function for the 
strange particle fraction of the hadronic gas, Z s , and we follow the notation of Ref. 
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However, given the recently measured values of the thermal parameters (which allow in 
particular for a rather high temperature T ~ 210 MeV) it is necessary to sum over some 
more strange hadronic particles than was done in Eq. ([/]) of Ref. [|l4j] . Including the possibility 
of only partially saturated strange phase space through the factor 7 S above, but neglecting 
for the moment the isospin asymmetry 8ii, we have: 



InZ 



HG 



VhT 3 
2tt 2 



(A s A q 1 + K%hsF K + (A S AJ + X^X^HFy 



+(A s 2 A q + A s - 2 A q 1 ) 7s 2 F H + (A s 3 + A s " 3 ) 7s 3 F fi 



(10) 



where the kaon (K), hyperon (Y), cascade (E) and omega (Q) degrees of freedom in the 
hadronic gas are included successively. The phase space factors of the various strange 
particles are: 



F K 


= J2g K W(m K jT) 

j 




= K,K*,K*,... 


, m < 1650 MeV , 


F y 


= J29y 3 W(m Y jT); 

j 


Yj = 


-- A, £,£(1385),. 


., m < 1750 MeV , 


F s 


= J29^W(m Bj /T); 
j 


= 


= 2,5(1530),..., 


m < 1820 MeV , 


Fn 


= J29n J W(m nj /T); 




= a 





where the are the spin-isospin degeneracy factors, W(x) = x 2 K 2 (x), and K 2 is the modified 
Bessel function. 

For the quark-gluon plasma, on the other hand, the strange contribution to the partition 
function is much simpler, because the strange quarks and anti-quarks are isolated and do 
not occur in clusters with non-strange quarks as in the hadron gas: 



In Z QGP = 9sVT^_ j p2 ^ ^ ^ + 7sV -V^V^ + ln ^ + 7sA - V V^/t 



with the strange quark spin-color degeneracy factor 



2x3 = 6. 



• (12) 



(13) 



For the value of the strange quark mass we take m s ~ 150-180 MeV; in Eq. (12) 7s and 
A s do not factorize as has been the case for the hadron gas, because T is of order m s and 
the Boltzmann approximation is not always satisfactory and can produce errors of the order 
of 20%. Recall that in general the strange particle density in a QGP is larger by a factor 
2-5 than in a HG, due to the lower threshold (m s < mx) and the presence of the color 
degeneracy factor in Eq. ([13|) . 



2.4 Strangeness balance 

We will now discuss the reduction of the number of free parameters caused by the condition 
of strangeness balance. Since the net strangeness of the colliding nuclei in the initial state is 



10 



J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank and J. Rafelski 



zero, and strangeness is conserved by strong interactions, the net strangeness of the collision 
fireball will stay close to zero throughout the collision, except for corrections due to possi- 
ble strangeness asymmetry in pre-equilibrium and asymmetric surface radiation of strange 



hadrons [22]. If these processes can be neglected, the resulting condition of vanishing total 



strangeness takes the form 



= p s = (s)-(s)=X s ^-\nZ s , (14) 



where p s is the net strangeness density, 

Ps = J2 S iPii ( 15 ) 

i 

with Si being the strangeness of particle species % having density p«, and the sum going over 
all particle species % existing in the respective phase (QGP or HG). While for the QGP 
this relation leads always to the result A s = 1 or /j s = (as one easily verifies by inserting 
Eq. fll2|)), for the hadron gas it is an implicit equation relating A s to A q in a way which 



depends on the temperature T [14], |23| . 

In order to also be able to account for pre-equilibrium and surface emission fluctuations, 
we introduce the net strangeness fraction 

£ = ^T' (16) 



and will also consider solutions for the analogue of Eq. ([14]) , where the left hand side is 
given by —s{s). In this work we restrict ourselves to asymmetries |e| < 0.1 which might be 
accessible by purely statistical fluctuations caused by surface emission, without "strangeness 
distillation" during QGP hadronization |f2"5| , |2"!2]j which arises when chemical equilibrium is 
enforced during a slow QGP hadronization process. An analysis of the phase diagram and 
the strangeness balance conditions for systems with very large net strangeness is presented 



elsewhere ||24|| . We note briefly the practical consequences of e = ±0.1: for the fireball 
created in central S-W collisions one has about 108 baryons in the interaction tube. If we 
consider a hadron gas with the thermal parameters (without flow component) as determined 
for these S-W collision data ||, we find a strange pair abundance of about 0.4 per baryon 
or about 40 strange pairs in total. The maximum asymmetry considered here in the HG 
phase at T ~ 200 MeV is thus about ±4 strange quarks. Fluctuations of this magnitude are 
expected to occur in about half of the collision systems. 



2.4.1 The case //„ = (A s = 1) 

Since the vanishing of p s for all values of T and /ib in a strangeness neutral QGP appears to 
be such a unique feature of this state, it is interesting to ask to what extent this condition 
could arise in a HG. We will therefore first study the mechanisms which lead in a HG model 
to a small strange quark chemical potential in a restricted region of temperatures. Of course, 
for a baryon-free system (/xg = 0) even in the hadron gas fi s = is a natural consequence 
of strangeness neutrality. For systems with non-vanishing baryon number (like the central 
rapidity fireballs created at AGS and SPS energies) this is no longer true. We therefore 
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ask ourselves under which conditions /x s — 0, viz. A s = 1, allows for non-zero values of /x B 
to be consistent with strangeness neutrality. This question can be answered analytically: 
inserting the partition function (^) into the strangeness neutrality condition (|l4|), we see 
that at A s = 1 the phase space for f2 and fl cancels out, giving the exact result 

"B = 3cosh-(^- 7 ,|). (17) 

There is a real solution only when the argument on the right hand side is larger than unity. 
It turns out that this condition is a sensitive function of the temperature T and of the 
hadronic resonances included (and therefore of the masses of the resonances). In particular, 
for any given resonance mass spectrum used to compute the phase space factors Fj, there is 
a temperature T beyond which no such solution is possible — this occurs since F^/Fy is a 
monotonically decreasing function of T (see Fig. 1 in Ref. ||T4||). 



The case of exact strangeness neutrality (e = 0) in a HG in absolute chemical equilibrium 
(7 S = 1) was considered in 0, and for this choice of parameters our following theoretical 
analysis essentially agrees with that work. 



2.4.2 HG behavior of /x s and /xb for nearly conserved strangeness 

An important aspect to consider is the value of /x s observed and how it is compatible with 
the picture of the reaction. The reason that /x s is of such importance is the expected value 
/i s = associated with QGP. We first will address the question under what circumstances 
this value is accessible to a HG fireball, and then turn to discuss more generally how different 
dynamical scenarios can be distinguished. We expect that there is a domain of temperatures 
for which even a considerable change of /xb between /xb = and /x B does not induce a 
significant change of /x s . 

For the extensive set of resonances included here this quite peculiar behavior occurs at 
T ~ 210 MeV for e = and at T ~ 200 MeV for e = -0.1. In Figs.ga, gb, |c we present the 
constraint between /x s and /Xb at fixed T = 200 MeV (solid curves), 150 MeV (long-dashed 
curves), 300 MeV (short-dashed curves) and 1000 (~ oo) MeV (dotted curves) for e = 
(Fig.|a), e = 0.1 (Fig. |b), and e = -0.1 (Fig.|c), all for 7s = 0.7. The choice % = 1 would 
not change our results significantly. The flat behavior of function of /xb for e = 

and T ~ 200 MeV makes it hard to distinguish the hadron gas from a QGP (which has 
/i s = independent of /xb) in this temperature range. This was pointed out in Ref. @, but 
we think that Eq. (|PT|) and the corresponding Fig.^|a, with our definition of fi s (see footnote 
|l|), exhibits this effect more clearly and allows us to understand the accidental nature of this 
behavior. 

For e = (Fig. |]a) all curves pass through /t s = at /ib = 0, and as T — > oo one ap- 
proaches a limiting curve which is rather well represented by the dotted curve corresponding 
to 1000 MeV. For e^Owe note that in the QGP phase the strange quark chemical potential 
does not vanish (see horizontal lines in Figures^b, 0c), and that there is a common crossing 
region of the HG and QGP results at a finite value of /x s , close to /xb = 0. For T — > oo for 
a fixed positive (negative) net strangeness fraction e, /x s continues to increase (decrease) at 
fixed /xb as a function of T, and only the slope as a function of /xb approaches a definite 
limit as T — > oo. 
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Figure 2: Strangeness chemical potential fl s ver- 
sus baryon chemical potential /ig . The long-dashed 
line corresponds to T = 150 MeV, the solid line to 
T = 200 MeV, and the dashed line to T = 300 
MeV. In Fig.|a, for e = (s - s) / S = 0, the dot- 
ted line is the limiting curve for large T and the 
horizontal solid line corresponds to /i s = in QGP 
independent of T. Fig.||b, c are the same as Fig. ||a 
with £ ± 0.1 except the dotted lines corresponding 
to T = 1000 MeV. The horizontal lines give /i s in 
QGP for the corresponding temperatures. 
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3 Analysis of the data within the fireball model 

In this section we describe, using as an example the experiment with the most detailed infor- 



mation of strange baryon and anti-baryon production, WA85 p6[ , how to extract estimates 
of the relevant thermodynamic parameters from particle production data. This procedure is 
straightforward if one can focus on the ground state baryons of a given flavor content, i.e. 
if the production and decay of heavier resonances can be neglected. Such an analysis was 
presented in || [J, and its essential steps will be shortly repeated in subsections 3.1 through 
3.3. If T is indeed in the region of 210 MeV as the m^-slopes seem to indicate, resonance 
production and decays can, however, not be neglected, and therefore such a procedure is not 



fully consistent. In subsection |3]4] we therefore redo this analysis including resonance decays, 
which is then much more involved. The main result is that resonance decays mainly affect 
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the extraction of the strangeness saturation factor 7 S (inclusion of resonance decays drives 
the extracted 7 S closer to 1, i.e. towards full saturation of hadronic gas phase space), while 
the striking conclusion that the data point towards a vanishing strange chemical potential 
/i s ~ is remarkably stable against these corrections. 



3.1 Cancellation of the Boltzmann factor 

It was argued in || that within the thermal model it is quite easy to extract the fugacities 
from a given particle ratio at fixed mj_ and y, thus eliminating the need to extrapolate the 
experimental data to full phase space. Specifically, for a stationary fireball one finds for 
the ratio of particles i and j, for each of which the invariant cross section is experimentally 
known at a certain rapidity y ± dy/2 in the region m±_ > mj^: 

9i li Xi ■& m -L cosh V exp(-m_L cosh y/T) dm\ 



dNi/dy 



dNj/dy 



>m c , 



9j lj Xj Im cut m J- cosh y exp(— m_|_ cosh y/T) dm\ 



Here we used the identity E = m± cosh y for the particle energy (with y measured relative 
to the center of mass of the collision fireball) to rewrite the transverse mass integral over the 
Boltzmann distribution E exp(—E/T) (the factor E results from starting with the Lorentz- 
invariant cross section E d 3 N/d 3 p). If the fireball is not stationary, but expands in the 
longitudinal direction in a boost-invariant way, dN/dy is independent of y, and the factors 
E exp(—E/T) under the integrals on the r. h. s. of Eq. flT8|) have to be replaced by a Bessel 
function |J: 

where we neglected a common normalization constant containing the fireball volume which 
drops out from the ratio!. The important observation is that obviously in both cases the 
Boltzmann factors (and with them all dependence on the particle rest masses) drop out 
from the ratio, due to the common lower cutoff on m_|_ (which, of course, must be larger 
than either rest mass). The same remains true if m± is not integrated over, but the ratio is 
evaluated at fixed m±. Since the degeneracy factors gi are known, this immediately yields 
the ratio of the generalized fugacities 7,Aj (see Eq. @). 



3.2 Extraction of fi s and /x q from anti-baryon data 

Thus, comparing the spectra of particles with those of their anti-particles within overlapping 
regions of m±, the Boltzmann and all other statistical factors cancel from their relative 



normalization, which is thus a function of the fugacities only. [In subsection 3^ we show 
that this statement, although no longer strictly true, still remains approximately valid when 
feed-down by resonance decays is included.] For the currently available two such ratios one 
has specifically 

i?H=="^, (20) 



2 The same result is obtained if the spectrum in Eq. ( |l8| ) is integrated over rapidity, i.e. if one wants to 
describe with a stationary fireball data obtained over a large rapidity window. 
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R A = — = ^±_K^ 1 (21) 
A AdA u A s 

These ratios can easily be related to each other, in a way which shows explicitly the 

respective isospin asymmetry factors and strangeness fugacity dependence. Eqs. ( p0| , pT|) 
imply: 

R A R^ = e 6 ^ T ■ e 25 ^ T , (22) 

R E Rf = e 6 ^ T ■ e-^' T . (23) 



Eqs. (^, |23|) are generally valid, irrespective of the state of the system (HG or QGP), 
as long as the momentum spectra of the radiated particles are thermal with a common 
temperature. We see that once the left hand side is known experimentally, it determines 
rather accurately the values of /i q , /i s which enter on the right hand side with a dominating 
factor 6, while the (small) isospin asymmetry Sfi plays only a minor role. This explains how, 



by applying these identities to the WA85 data it has been possible to determine the 



chemical potentials with considerable precision in spite of the still pretty large experimental 
errors on the measured values of R A , R^. 

Turning to the pertinent experimental results we first recall the recently released data 
on the rapidity dependence of A and A production rates measured by the experiment NA36 



27]. There is a clear indication that the production of these particles in nuclear S-Pb 
collisions is concentrated in the central rapidity domain — it does not follow in its behavior 
the FRITIOF simulations which, on the other hand, are consistent with the production rates 
from p-Pb collisions. This supports the notion of a central fireball reaction picture for S-W 



collisions studied by WA85 |[26|| . The most striking result of the latter experiment is the fact 
that the abundance of charged anti-cascades (ssd) is unusually enhanced, in particular when 
compared to the abundance of cascades (ssd) and anti-lambdas (sud). We now interpret the 
data in terms of the local equilibrium model, thus fixing the values of /x s and /i q . 
The ratio has been reported as: 

R s = 0.39 ± 0.07 for y G (2.3, 3.0) and m± > 1.72 GeV. (24) 

Note that, in p-W reactions in the same (p±,y) region, a smaller value for the i?= ratio, 
namely 0.27 ± 0.06, is found. The A/A ratio is: 

R A = 0.13 ± 0.03 for y G (2.3, 3.0) and m± > 1.72 GeV. (25) 

In Eq. corrections were applied to eliminate hyperons originating from cascade decays, 
but not those originating from decays ofQ—>-A + KorQ—>-A + K which are of little 
significance for the high m± considered here. The ratio R\ for S-W collisions is slightly 
smaller than for p-W collisions in the same kinematic range. 



From these two results, together with Eqs. (^, |23|, and the value Rf = 1.08 mentioned 
after Eq. (0), we obtain the following values of the chemical potentials for S-W central 
collisions at 200 GeV A: 



(26) 
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^=f^(Rf-l) = 0.042 ± 0.008 , (27) 

/i s hiRJRl -0.084 . , 

t£ = = -0.018 ±0.05. 28 

1 6 

A study of the temperature observed in these collisions was recently made ||. As noted 
in || and discussed below in section 4, the possibility of a HG interpretation of these data 
only poses itself if the transverse mass slope of the produced particles is entirely due to the 
thermal motion. Under this assumption the temperature is T = 210 ±10 MeV, and therefore 
we have: 

// B = 335 ± 25 MeV , 6/jl = 9 ± 2 MeV , /t s = -3.8 ± 10 MeV . (29) 

The last result translates into the value A s = 0.98±0.05 for the strange quark fugacity. Much 
of what we do later in this paper will rest on this stunning result that the strange chemical 
potential is very small and perfectly compatible with zero. In our calculations below we 
will thus consider in particular a HG constrained to this range of values (as well as to near 
strangeness neutrality) . 



3.3 Determination of 7 S 



A complete cancellation of the fugacity factors occurs when we form the product of the 
abundances of baryons and anti-baryons. We further can take advantage of the cancellation 
of the Boltzmann factors by comparing this product for two different particle kinds, e.g. we 
consider: 



A A 



(30) 



If the phase space of strangeness, like that of the light flavors, were fully saturated, the fireball 
model would imply T s = 1. However, any deviation from absolute chemical equilibrium as 
expressed by the factor 7 S will change the value of T s . Ignoring as before any feed-down 
effects from higher resonances we find 



7 S 



(31) 



In principle the measurement of 7 S can be done with other particle ratios; in the absence of 
resonance feed- down we have 



7 S 2 



A A 

p p 



m±>m c l 



A A 



mj_>m c . 



22- 22- 



(32) 



m±>m c . 



where in the last relation the factors 2 in the denominator correct for the spin-3/2 nature 
of the fl. As we will see in subsection 3.5, Eq. (0), these double ratios are considerably 
affected by resonance decays. 

We note that in the kinematic domain of Eqs. ( |24"| , |25|) the experimental results reported 
by the WA85 collaboration are: 



(33) 



16 



J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank and J. Rafelski 



If the mass difference between A and S° is neglected, this implies in the framework of the 
thermal model that an equal number of A's and S°'s are produced, such that 

^ = 1.2 ±0.4, =- = 0.40 ± 0.08 . (34) 

A A 

The fact that the more massive and stranger anti-cascade exceeds at fixed m± the abundance 
of the anti-lambda is most striking. These results are inexplicable in terms of simple cascade 
models for the heavy-ion collision ||28|| . The relative yield of H~ appears to be 5 times greater 



than seen in the p-p ISR experiment and all other values reported in the literature p6 |. 
Combining the experimental result Eq. ( |3~4"1) with Eqs. (^, |3l|), we find the value 



7 S = 0.7 ±0.1. (35) 
3.4 Temperature, transverse mass spectrum 

The data from the WA85 collaboration |2(| correspond to an apparent temperature (from 
the slope of the m_|_-spectra) of T a p parent = 210 ± 10 MeV ||. In our numerical studies we 
therefore concentrated on two extreme cases: 

i) T = 210 MeV, (3±_ = (no transverse flow), which we call the "thermal" scenario, 

ii) T = 150 MeV, /3± = 0.32, which we call the "flow" scenario for which the freeze-out 
temperature has been estimated following the results of M. 

If the slope of the experimental spectra is not entirely due to thermal motion, but con- 
tains a flow component due to collective expansion of the emitting source |J, then "flow 
spectra" (i.e. boosted Boltzmann spectra) have to be used in the numerator and denomina- 
tor of Eq. (|i~8P. If the flow is azimuthally symmetric and boost-invariant in the longitudinal 



direction, the thermal spectrum (H5) has to be replaced by (again neglecting a common 
normalization constant) 

oWf° w . „ (m ± coshp\ (p ± smhp\ 

d^ ± = 9i 7 * Xl m± Kl (— r— J /o I" T - J ' (36) 

where K\ and Iq are the modified Bessel functions and 

p = tanh -1 (3 f (37) 

is the transverse flow rapidity. This spectrum has an asymptotic slope corresponding to an 
apparent "blue-shifted" temperature 



-^apparent T A 



The same result as in Eq. fl3"6] ) is obtained if the spectrum from a spherically expanding 
fireball is integrated over rapidity. Since some longitudinal flow is likely to be present in 
S-W collisions, and also the kinematic window of the WA85 experiment is nearly one unit 
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of rapidity, we think that the form Q35|), in which integration over a large rapidity window is 
implied, is most appropriate for an analysis of the experimental spectra. (The "large window" 
combines the one unit of ^-acceptance in WA85 with about 1.5 units of flow expected || for 
this system.) 

Obviously, through the separate dependence on m± and p±_ Eq. ( |36|) depends on the 
particle rest mass explicitly. This spoils the cancellation of the spectral shape from the 
m_i_-integrated particle ratios as well as their independence from the m_i_-cut, and thus also 
renders the extraction of the fugacities, etc. more complicated. 

In order to assess the importance of these effects we devote next subsection to a more 
complex reanalysis of the experimentally measured strange baryon and anti-baryon ratios 
which takes all these complications into account. Anticipating this development we show 
in Fig. |^ the m^-spectra of the various particle species, for both the "thermal" scenario 
(Figs. 0a, b) and the "flow" scenario (Figs.|3]c, d). Dashed lines indicate the direct thermal 
contribution to the spectra, while the solid lines include all resonance decay contributions 
(see following subsection). These results allow to directly assess the importance of resonance 
decays on the m^-slopes of the observed particles and on the particle ratios. For the Q there 
are no resonance contributions, due to the imposed mass cutoff in the resonance spectrum. 
While the overall normalization of the spectra is arbitrary, their relative normalization is 
according to the fugacities given in Table 1. 

The presence of flow affects the spectra mostly at low m± near the mass threshold 
[[J, where it leads to a flattening, i.e. to a higher apparent temperature. The effect is 
proportional to the particle mass (since for small p± the effect of flow can be estimated by 
p flow ~ m fflow) and thus most clearly seen in the S~ and Q spectra. 



3.5 Resonance decays 



As already mentioned, Eq. ([18]) is only correct if feed-down to the observed particles by 
strong decays of higher resonances is neglected. At temperatures of order 200 MeV these 
resonance decays can no longer be totally ignored. This has two annoying consequences: 
first, since the decay products have a steeper m^-spectrum than the original resonances Jl 



the various contributions to numerator and denominator in ( |18|) have different slopes, and 
the integrals over the Boltzmann factors no longer cancel exactly. Second, since different 
sets of resonances contribute to i and j, the ratio is modified by an a priori unknown factor 
which turns out to be hard to calculate and in general does depend on the experimental 
m^-cutoff as well as on the temperature. 

Since the decay spectra drop more steeply as a function of m± than the thermal ones, 
the role of resonance feed-down can in principle be suppressed by going to sufficiently large 
m±. For pions and kaons, for example, the thermal contribution dominates for m± > 1-1.5 
GeV, in spite of the huge contribution of decay pions to the total yield |HJ; the latter all 



end up at lower m±. For baryons due to their larger masses, the difference in slope between 
the thermal and decay baryons is less drastic, and one has to go to still larger m± before the 
resonance decay contributions are sufficiently suppressed. 
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Including resonance decays, the numerator and denominator of Eq. (ITS) take the form 



dy 



r dm 2 f rfAf ermal/flow (T) _ 6 dN?(T) \ 
Jmf* ~ L 1 dy dm\ R ~* 1 dy dm\ J ' 



with the direct thermal or flow contribution given by ([19]) or fl36|), respectively, and the 
contribution from decays of resonances R —> i + 2 + . . . +n (with branching ratio into 
the observed channel i) calculated according to [F7| 

AA[R rs+ r Y+ r M+ 

' dsg n (s) / dY / _ dM\ (40) 



dydm]_ J S - Jy- Jm 

M ( dN R 



P\p\ - [ME* - M ± m ± cosh(F - y)] 2 \ dY dM ±J 

Capital letters indicate variables associated with the resonance R, y 7 ? is the invariant mass 
of the unobserved decay products 2, . . . , n, and the kinematic limits are given by 

s Jj2m k ) , S+ = (M-m 4 ) 2 , (41) 



\k=2 



Y±=y ± sinh" 1 — , (42) 
m± 



_E*m± cosh(y — y) ± p±\Jp* 2 — m 2 ± sinh 2 (y — y) 
m 2 L sinh 2 (y — y) + m 2 



Mt=u ^.^^,^ „ ^ (43) 



E "=Wi {M ' + ™' " s) • p ' = ^ E " 2 - ,n >- < 44) 

In (|40D <7 n (s) is the decay phase space for the n-body decay i? — > i + 2 + . . . + n, and for the 
dominant 2-body decay (assuming isotropic decay in the resonance rest frame) it is given by 

g 2 (s) = J-5(s-M 2 ), (45) 
Anp* 

for n > 2 see Ref . • 

The resonance spectrum dN^/dY dM\ entering the r. h. s. of Eq. ( f4"0"D is in general 
itself a sum of a thermal or flow spectrum like ( fL9|) or (^) and of decay spectra from still 
higher lying resonances; for example, most decays of high-lying nucleon resonances into 
protons or anti-protons proceed through the A(1232) or its anti-particle (e.g. iV(1440) — > 
A(1232) + 7r — ► p + 2ir). This has been taken into account in our calculations. Extracting 
the degeneracy factors and fugacities of the decaying resonances, we write shortly 

poo dN^ 

A* s 7R X R N» . 7R A H / *»» ? 3R ^—^ , (46) 



with the sum now including only resonances with identical quantum numbers; if these quan- 
tum numbers agree with those of species i, the sum is meant to also include the thermal 
contribution. 
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Between particles and anti-particles we have the relation 

Nf = ~ fR \- R 1 N l R = \- R 2 N l R . 



The three independent particle ratios measured by WA85 [26[] are thus given by 



A 
A 



m i >m t ; ut 



%VA^vr + 7.yA g - 2 jyf + 7 S 3 A s - 3 A>r 

7s A^A s iVl* + 7s 2 A q A s 2 iVf * + TfA^iVf 

7 b 2 A q A b 2 iVr + TjA|iVr 
7sA 2 A s A>I* + 7s 2 A q A s 2 iVf * + 7s 3 A s 3 iVf ' 



(47) 

(48) 
(49) 
(50) 



N\* contains also (in fact as its most important contribution) the electromagnetic decay 
S° — > A + 7, and N^* includes the weak decay Q —>■ YK (since fl's have not been recon- 
structed so far by WA85). If the right hand sides are evaluated using Eqs. ([50, |4(], ffl) with 
m cut = x.72 GeV, for the left hand sides the experimental values (|24 
and Eqs. 



31) can be inserted, 



|, |49|, |5(J) can be solved for 7, 
between A u and Ad since this would require the calculation of a separate decay cascade for 



s , A s , and A q . Note that we do not distinguish here 



each member of an isospin multiplet. The smallness of 8fi in (27) justifies this approximation. 

The results of our analysis are given in Tables 1 and 2 and Figures || and ^. Comparing 
the results for the thermal parameters in Table 1 with Eqs. 



28| , j35| ) which contained only 



an estimate of the radiative E° decay but no contributions from higher lying resonances, we 
can draw the following conclusions: 



1. The extraction of A q and A s according to Eqs. (20, 21) is only very weakly affected by 
resonance decays and by the origin of the slope of the m^-spectrum (thermal or flow). 
(The absolute value of the associated chemical potentials does, of course, depend on 
the choice of the freeze-out temperature.) The reason is that in Eqs. ( |48"D and ( f49|) 
in each case the first term in the sums occurring in the numerator and denominator 
completely dominates. (For S~'s the effects from fi-decays are small, and for A's 
the strong S° contribution overwhelms the effect from 5* decays.) Taking the results 
for the "thermal" scenario from Table 1 and evaluating with them the A/A, 
and S~/A ratios excluding the contributions from higher resonances changes them to 
0.383, 0.121, and 0.262, respectively (instead of the input values 0.39, 0.13, and 0.20, 
respectively). Leaving out also the S° — > A7 contribution changes only the 5~/A 
ratio which now becomes 0.430. Obviously this ratio is the one which is most strongly 
affected by the resonance decays. 

2. Since the H~/A ratio is the crucial ingredient for the determination of the strangeness 
saturation factor 7 S , the effects from resonance decays and flow can be clearly seen in 
this number. From Eq. (0) we had extracted 7 S = 0.7; improving the rough factor 2 
estimate for the S° contribution by its correct thermal weight at T = 210 MeV would 
have lowered 7 S to 0.57. The decay of higher- lying resonances raises the value of 7 S up 
to 0.76 at T = 210 MeV. If due the presence of flow the m_i_-slope does not give the true 
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(a) pure "thermal" 
T = 210 MeV, P f = 


(b) "thermal & flow" 
T = 150 MeV, /3/ = 0.32 


As 

/i s (MeV J 


0.973 
-0.027 

— 5. ( 


0.972 
-0.028 
— 4.Z 


K 

/i q (MeV) 


1.718 
0.54 
113.6 


1.710 
0.54 
80.5 


7s 


0.756 


0.902 



Table 1: Thermal fireball parameters extracted from the WA85 data |26) on strange baryon and anti-baryon 
production, for two different interpretations of the measured m±-s\ope. Resonance decays were included. 
For details see text. 



1 J ■m ± >m c ^ t 


(a) "thermal" 


(a) "thermal" 


(b) "flow" 


(a) "flow" 




m c f = 1.72 GeV 


m]_ ut = 1.2 GeV 


ml ut = 1.72 GeV 


m]_ ut = 1.2 GeV 




0.390* 




0.390* 




A/A 


0.130* 


0.135 


0.130* 


0.133 


S-/A 


0.200* 




0.200* 




S-/A 


0.600* 




0.600* 




A/p 


0.625 


0.617 


0.888 


0.867 


A/p 


2.044 


2.084 


2.833 


2.802 


n~/Er 


0.498 




0.311 






1.502 




0.944 




p/p 


0.040 


0.040 


0.041 


0.041 


K°/A 


0.086 


0.077 


0.138 


0.152 


K°Jp 


0.054 


0.048 


0.123 


0.132 



Table 2: Input values (*) and predicted particle ratios for the WA85 experiment. The parameters for the 
"thermal" and "flow" scenarios are as in Table 1. 
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qf the spectra is arbitrary, their^kitive 
lines indicate direct thermal prodis 
Qar^ge decays. For parameters and defaSi; 



temperature, tBS? estimator 7 S go<f§°up even higher: WS T = 15eP<MeV, (3^0.32 (i.e. 

m 1 [MeV] m 1 [MeV] 

values which are consistent with the kinetic freeze-out criterion developed in |J) we 
find 7 S = 0.9, c.f. Table 1, corresponding to nearly complete saturation of the strange 
phase space! 

3. The conclusion that the WA85 data on strange baryon and anti-baryon production 
in 200 GeVA S-W collisions seem to indicate a vanishing strange quark chemical 
potential is extremely stable with respect to resonance decay contributions and the 
possibility of a flow component in the spectra. It will serve as the cornerstone of the 
following studies in section 4. 
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Figure 4: Various 77l_i_-integrated particle ratios as a function of the low-77lj"*. Solid lines are for the 
"thermal" scenario, dashed lines are for the "flow" scenario, with parameters as given in Table 1. 



3.6 Particle ratios 

It is clear that the extraction of three thermodynamic parameters (A q , A s , and 7 S ) from three 
independent particle ratios cannot be considered a very convincing test of the thermal pic- 
ture, in particular since the extracted values depend on the choice of the temperature whose 
independent derivation from the shape of the m^-spectrum is not quite unique (unless one 
accepts the theoretical consistency arguments of Ref. |9| which prefer the lower temperature 
of around 150 MeV). It is therefore mandatory that the picture receives further experimental 
tests in the form of other independent particle ratios. 
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In Table 2 and Figure |3] we have therefore included predictions for several additional 
particle ratios, in particular involving Q, Q and p, p which are expected to become soon 



available from the WA94 collaboration fl3(J in a kinematic region compatible with the WA85 
results. Our calculations explore the consequences of both the "thermal" and "flow" scenarios 
above, and anticipate also a possible improvement of the experimental acceptance towards 
lower transverse momenta, by calculating in Table 2 some of the ratios also with an m± cut 
of only 1.2 GeV, or plotting them (in Fig.|j) as a function of this cut. 

Comparing mj_-integrated particle ratios with different m± cuts gives an idea of the effects 
coming from the somewhat different slopes of thermally emitted particles and resonance 
decay products. From Fig. [| we see that these effects are generally small for ratios involving 
only baryons: due to their heavy mass, the slope difference of parent and daughter particles 
in resonance decays is small |17|, H, resulting in a slight steepening of the spectra at low and 
intermediate m±. For kaons, on the other hand, the decay contributions are considerably 
steeper than the thermal one, leading to an appreciable upward curvature of the spectrum 
at low mj_, while flattening out to the asymptotic apparent temperature for m± > 1-1.5 
GeV, where the thermal component dominates. In consequence, as seen in Fig.|4], the K/A 
and K/p ratios are proportionally more sensitive to the m± cut. 

It should be kept in mind that the large QfE~ and (in particular) ratios are partly 

due to the additional spin degeneracy of the spin-3/2 Q baryon. The apparent strong rrf^- 
dependence of all ratios involving Q or f2 baryons (especially in the flow scenario) should be 
taken with some caution: it results mostly from the absence of resonance decay contributions 
to the Q spectrum in our calculation and might be modified if (so far unknown) Q* resonances 
existed and should be included. 

When comparing different ratios involving A's (e.g. in order to estimate the systematic of 
these ratios with increasing strangeness) , one should always keep in mind the considerable S° 
contribution which, at T = 210 MeV, effectively changes the degeneracy factor associated 
with the A by about a factor 1.6. Once this is corrected for, i.e. the 2° contribution is 
accounted for, one finds the inequalities 

2 AA 1 nil 

% > — > —T=r- > 7 =-S= ■ (51) 

pp AA 4 S-E- 

Paradoxically, Eq. ([H]) suggests an increasing importance of resonance decay corrections with 
increasing strangeness content, as compared to the estimate, Eq. ([32]) . This is due to the 
fact that the stranger baryons are also heavier and thus receive less resonance contributions 
(because the corresponding excited states are strongly suppressed by the mass term in the 
Boltzmann factor) than the lighter ones (whose resonances are easier to excite). In part this 
tendency may be somewhat exaggerated by our cut in the resonance spectrum at about 2 
GeV; taking into account even higher resonances (although it is not clear whether can indeed 
be excited with statistical probabilities in nuclear collisions) might modify in particular the 
fi/E - ratios. 

We regard the ratiosi predicted in Table 2 and Fig. 01 as a qualitative test of the idea 



3 Ratios involving protons should be taken with a grain of salt since it may not be easy experimentally 
to remove the contamination by cold projectile and target spectator protons. Similarly, those ratios which 
involve kaons may come out differently than predicted here because, as we argue in section the assumption 
of chemical equilibrium may break down for kaons if they are created by hadronization of a QGP. 
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that the emitter is a fireball which is in thermodynamic equilibrium with the given thermo- 
dynamic parameters. However, even if this test is eventually passed, the fact that we can 
successfully describe all observed particle ratios in this way does not yet tell us immediately 
the nature of the emitting source (HG or QGP). That information resides in supplementary 
information concerning the global properties of the source and their relationship to the ex- 
tracted thermodynamic parameters, as well as in the systematic behavior of these parameters 
under varying experimental conditions (collision energy, projectile and target size, collision 
centrality etc.). These important issues will be discussed in more detail in the following 
sections f| and §. 

3.6.1 Relation between strange baryon ratios from the constrained fireball 

Since thermal particle ratios are functions of only fis and fi s , the accidental vanishing of /j, s 
at a temperature T ~ 200-220 MeV due to the constraint to strangeness neutrality removes 
the opportunity to simply distinguish the two phases — a HG phase of strongly interacting 
matter, at this level of discussion, is indistinguishable from QGP as long as a measurement 
of the relative multiplicities of strange particles is considered at one, and only one collision 
energy. Even in this case many theoretical questions arise, such as why is /i s = or why is 
7 S near to unity? However, as Figs. 0a, ^|b, ^c clearly show, this will not remain a problem at 
higher or lower temperature. To make this point more quantitative, we now study strange 
baryon ratios arising at different temperatures from thermal fireballs constrained to a (nearly) 
vanishing strangeness. 

In Fig. [| we show, for the case of exactly vanishing strangeness and for various temper- 
atures, the resulting relation between R^, and with R\. In addition to the HG results for 
temperatures T = 200 MeV (solid line), T = 150 MeV (dashed line) and T = 300 MeV 
(dotted line) we show the case /z s = corresponding to a QGP source hadronizing rapidly 
by quark recombination (dashed-dotted line). The cross corresponds to the result reported 
by the WA85 experiment ||26|| . As can be seen, the /x s = curve nearly coincides with the 
T = 210 MeV curve, as noted before |], 0j, but fireballs decoupling at significantly lower or 
higher temperatures (as expected at lower and higher collision energies, respectively) would 
lead to significantly different particle ratios if the thermal model is correct. 

3.6.2 High-m^ kaon abundances 

In addition to the baryon ratios we can also consider the ratio of kaons to hyperons, again at 
fixed m±. Because of the experimental procedures used, which rely on the observation of the 
disintegration of neutral strange particles into two charged decay products, a comparison of 
the K® with the A (which includes the A's from the electromagnetic decay of E°'s) is the 
most useful next step. We therefore introduce: 





(52) 



A + S° 



where the second identity is again, of course, only valid if resonance decay contributions can 
be neglected. For the rather light kaons there exist, however, many different possibilities 
for secondary production through resonance decays, which limits somewhat the practical 
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Figure 5: Ra versus The long-dashed line 

corresponds to T = 150 MeV, the solid line to 
T = 200 MeV, and the dashed line to T = 300 
MeV in the HG. The dashed-dotted line corre- 
sponds to QGP (yU s = 0). 



Figure 6: R\ versus Rk with the same conven- 
tions as in Fig.||. 



usefulness of this ratio. For example, if we use the thermodynamic parameters extracted 
in subsection [3.4| from the WA85 data, we obtain for the "thermal" interpretation of the 
m^-slope (T = 210 MeV, no flow, m± = 1.72 GeV) a value Rk = 0.086, if all resonance 
contributions are included, while the ground state contributions alone according to Eq. fl52]) 
would give Rk = 0.102 (Rk = 0.125 if the A-£° mass difference is taken into account). 
Since the size of these corrections depends on the fireball temperature and the mj_-cut, a 
full quantitative study does presently not appear practical. 

Given these uncertainties, and in order to demonstrate the insensitivity of the iT/A-ratio 
to the nature of the fireball, we show in Fig.|6] (as a matter of principle) the uncorrected ratio 
Rk as a function of the uncorrected ratio R\: as the figure clearly shows, this observable is 
rather insensitive to both the temperature of the source and to its intrinsic structure (/i s = 
or fi s ^ 0). The lines are as in Fig. |5], and we note again the near coincidence of the T = 200 
MeV HG result with the QGP result. But even as T is changed to 150 or 300 MeV, the 
variation of R\ at fixed i?K stays well within the present experimental error bars and will 
thus not permit to identify the properties of the source. 



4 Total particle multiplicity and entropy 

We now return to discuss in detail a very important observation which is related to the 
particle multiplicity and charge flow associated with strange baryon production | pTj| , |5|. Even 
though no data are yet available on multiplicity densities in coincidence with strange particle 
production, the properties of the HG and QGP fireballs with regard to their entropy content 
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are so drastically different that we can draw at least qualitative conclusions by comparing 
multiplicity data from S-Pb collisions with S-W strangeness production, despite the fact 
that the two targets ( 207 Pb and 184 W) differ slightly in mass. 

QGP and HG states are easily distinguishable in the regime of values /ia, T, % of interest 
here. The entropy per baryon in the HG is iS HG /£> = 21.5 ± 1.5. Consequently, the pion 
multiplicity which can be expected from such a HG fireball is at best dN/dy = 4 ±0.5. This 
is less than half of the QGP based expectation, which we can easily estimate by considering 
the perturbative QGP equation of state: Up to corrections of order (fiB/nT) 2 ~ 0.03, the 
leading term of the specific entropy from light quarks and gluons is: 

s qgp 3vr 2 ,T\ r 14(1 — 50a s /2l7r) 32(1 - 15a s /4vr) 



B 2 V/VL 15(l-2a s /7r) 45(1 - 2a s /vr) 

= 17 (T/fJL q ) for a s = 0.6 . (53) 

Adding to this the considerable entropy content of the strange quarks (neglecting for technical 
reasons the a s -corrections in this case), we find for the QGP a specific entropy of about 
45 units at the same value of T/fi q . Poor knowledge of the value of a s and the strange 
quark mass in the QGP phase make this estimate uncertain by about ±15%. Still, the 
difference to the HG result is considerable in terms of experimental sensitivity. Checking the 
theoretical sensitivity we note that the point at which the entropy of HG and QGP coincide 
and strangeness vanishes and A s ~ 1 is at T ~ 135 MeV, /xb — 950 MeV, quite different 
from the region of interest here. We therefore consider now an observable which is capable 
to distinguish between the HG and QGP via their entropy content. 



4.1 EMU05 data 

Total charged particle multiplicities (excluding target/projectile fragments) above 600 in the 
central region have been observed by EMU05 fTl] in 200 GeV A S-Pb collisions, correspond- 
ing possibly to a total particle multiplicity of as large as 1000. We believe that such large 
multiplicities are implied by a QGP scenario for the central fireball, while being hardly com- 
patible with the corresponding HG interpretation. The data of the EMU05 collaboration 
[Ll| were obtained with a high multiplicity "trigger" which makes them useful in comparison 
with the WA85 results which are also triggered on high multiplicity. Some of these charged 
multiplicity data from emulsions are shown as a function of rapidity in Fig.[7|. We depict 
the quantity Dq, the difference in the number of positively and negatively charged particles 



normalized by their sum [11 



D = !1 1I_ (54) 

The data points correspond to 15 fully scanned "central" events of 200 GeV A S-Pb 
interactions, with the trigger requirement being a total charged multiplicity > 300. Reaction 
spectators (target fragments) are not observed in this experiment. At central rapidity a 
value of 0.08-0.09 is found for Dq. We note the rise of Dq in the projectile and target 
pseudorapidity regions, indicating a more rapid decrease of the particle multiplicity (mostly 
pions) than of the positive particle excess (protons). Obviously this feature proves the partial 
transparency at CERN energies of the Pb nucleus to the incoming sulphur projectile; it will 
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be quite interesting to see if this feature persists for heavier projectiles and in particular how 
the shape is in Pb-Pb collisions. 



28 



J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank and J. Rafelski 



4 6 
pseudorapidity 



/ <- \, = 0.95 




/tij [MeV] 



Figure 7: Emulsion data for charged particle mul- 
tiplicity as function of pseudo-rapidity: the differ- 
ence of positively and negatively charged particles 
normalized by the sum of both polarities. (Curtesy 
of EMU05 collaboration, Y. Takahashi et al. 0). 



Figure 9: Entropy per baryon S/B as a func- 
tion of -Dq 1 for fixed A s = 1 and conserved zero 
strangeness. The lower, dashed curve does not take 
into account the strange hadrons. 



Figure 8: Dq as a function of for fixed 
A s = 1 i 0.05 and conserved zero strangeness in 
HG. 




In our theoretical work we address primarily the central region, in which the smallest 
values of Dq indicate the highest multiplicity production. 



4.2 Significance of Dq 

The numerical evaluation of Dq in a HG at A s = 1 ± 0.05, where we fix the temperature for 
each fis such that strangeness is conserved, is shown in Fig.||. We see that approximately 



Dr 



1.3 GeV 



for /i B < 0.6 GeV. 



(55) 



This HG result is extremely simple, considering the complexity of the calculation. We thus 
see that in the only permissible HG scenario for the strangeness source which (according 
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to Eq. has to have ^ ~ 0.34 GeV, Eq. (^) would predict a value of Dq ~ 0.26; 

this is incompatible with the EMU05 results, see Fig. [7], where the central value is Dq(y = 
2.5 ± 0.5) = 0.088 ± 0.007. At this value of D Q the HG phase would have /x B = 115 MeV, 
rather than the 340 MeV we have extracted from the strange baryon data. Despite the fact 
that both measurements use different targets we believe that we can draw the conclusion 
that a HG fireball is incompatible with the combined EMU05 and WA85 data. 

Physically the variable Dq represents essentially the inverse of the specific entropy. The 
measured very small value of Dq thus implies a very large value for S/B. In Fig.|9| we show 
the specific entropy as a function of Dq 1 for a zero-strangeness hadron gas at A s = 1. We 
see the very strong direct correlation between these two quantities. We supplement this 
observation in Figs.pUla, |10"lb, and [10| c by also showing S/B as a function of Dq 1 at fixed 



T = 150, 200, 300 MeV, varying e = 0, ±0.1. We observe that, within the range of these 
parameters, the small observed value of Dq always requires a rather large specific entropy 
S/B = 50 ± 4. Thus Dq 1 effectively measures the entropy per baryon, independent of the 
specific thermal conditions of the source. 

Please note that in this connection it is essential to correctly account for the strange 



particles. To illustrate this point and to understand the value of -Dq we show in Fig. |TT] the 
product Dq-(S/B) as a function of /ie for fixed A s = 1 ± 0.05, where akin to Fig.|9|, the 
lower set of curves was obtained without strange particles. We see that strange particles 
have considerable impact on the result. 

Both the limiting value and near constancy of the result obtained without strange par- 
ticles can be easily understood since in this case Dq-(S/B) is effectively the entropy per 
pion. To see this note that in absence of strange particles we have (neglecting the small u-al 
asymmetry and assuming pion symmetry N n + = N n - = N^o = N n /3): 

D « - °' 75 f TTTlkwK (56) 

where the last term in the denominator involves the sum over all initial charged particles 
other than pions, in particular heavy mesons, protons and anti-protons, etc. Considering the 
product of Dq with S/B, the baryon content cancels and the result is effectively entropy per 
pion renormalized (as it turns out numerically with a factor 2) by the degeneracy factors and 
resonance contributions. This also explains why the product is rather constant as fi^ (and 
along with it T) changes: there is a slow increase in the entropy per pion with increasing 
/iB since the resulting decrease in T enhances the relative importance of the heavier baryon 
resonances with the associated larger entropy per particle. 

The question arises whether the QGP fireball model would be in better agreement with 
the combined WA85 and EMU05 data, or if the disagreement of the HG picture with the 
data is rather an indication for the failure of the fireball model in general. In order to be 
able to answer this question in detail we need in principle to have a sophisticated model for 
the hadronization of a QGP. However, as long the final state is one with thermal momentum 
distributions, the relationship between Dq and S/B is likely to be similar to that of Fig. |9| 
even if the observed hadrons result from a hadronizing QGP. The reason is that the product 
Dq • {S/B) (in other words, the entropy per pion) will not deviate drastically from its equi- 
librium hadron gas value even in a state which is considerably out of chemical equilibrium. 
Noting that there is little entropy production in the QGP hadronization process p| we can 
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Figure 10: a) Entropy per baryon as a function 
of Dq 1 at fixed temperature (T = 150, 200, 300 
MeV) for conserved zero strangeness (same conven- 
tions as in Fig. Ha), b) and c) Same as Fig.[ic|a but 
for e i 0.1 (compare to Fig.|]b, c). 
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thus subsume that after hadronization the result would be a point close to the curve shown 
in Fig. p. We see that the value Dq ~ 0.09 implies S/B ~ 50, which is indeed of the order 
required by a hadronizing QGP with the thermal parameters extracted in section 3. We note 
that entropy per baryon is a dimensionless function of the dimensionless quantity T/fi-g and 
hence different QGP scenarios with the same T/fi-Q all have the same value of the entropy 
per baryon. 

As the result of this study we conclude that the charge flow ratio with the charged particle 
multiplicity is an excellent experimental measure of the entropy per baryon in the source, 
and that the current results are suggestive of a much more entropy-rich emitter at central 
rapidity than can be expected in hadronic gas models. 
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4.3 NA35 data 

Even though the use of thermal model is 
less convincing for the small collision sys- 
tem S-S, we will shortly consider the re- 
sults obtained by the NA35 collaboration 
34], H in S-S interactions at 200 GeVA. In 
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Figure 11: The product Dq ■ {S/B) as function 
of /iB for fixed A s = 1 ± 0.05 and zero strangeness 
(compare to Fig.||). The lower set of curves was 
obtained without including strange particles. 



this case we can base our discussion on the 
reported charged particle rapidity densities 
pjfl and the measured A/A ratio f34|. S-S 
interactions show at central rapidity visibly 
less baryon number stopping than S-W col- 
lisions; from the data on d(N + — N~)/dy 
and dN~ jdy in |35| we read off a central ra- 
pidity value of £>Q +S (y = 3) = 0.065. The 
raw A/A ratio (uncorrected for S-decays) 
is strongly peaked at central rapidity and 
rises there to a value A/ A = 0.34 (with er- 
ror bars of order 50%) [[34]]. The observed 
transverse mass spectra can again be inter- 
preted in terms of a thermal, longitudinally 

expanding central fireball with T ~ 200-220 MeV |L7| or, including a transverse flow com- 
ponent, by T ~ 150 MeV combined with a flow velocity (3{ ~ 0.3 B. 

Following Eq. ([H]), the asymmetric A/A ratio implies a non-zero combination of chemical 
potentials, /iB + (3/2) /x s ~ (0.8±0.4) T. We note that the experimental ratio decreases when 
averaged over a larger rapidity interval, implying even larger chemical potentials. If T ~ 210 
MeV, both the HG and QGP model imply /i s ~ for strangeness neutral systems, and we 
thus obtain /iB — 170 ± 85 MeV. For the "flow" scenario with T = 150 MeV, strangeness 
conservation in a HG allows the above relation to be satisfied with /x s ~ 17 MeV at fi-g ~ 95 
MeV, again with rather large experimental errors. 

For the "thermal" case we can use Fig. || to translate /iB into Dq; with /iB = 170 ± 85 
MeV we obtain Dq ~ 0.13 ± 0.06. Although the central value of Dq extracted in this way is 
twice the measured value for this ratio, due to the large experimental uncertainty of the A/A 
ratio the two numbers still agree within one standard deviation. We note that a HG with 
Dq = 0.13 possesses a specific entropy S/B = 35, while Dq = 0.065 requires S/B = 65-75. 
For the "flow" scenario, on the other hand, we compute with the conventional HG-EoS (with 
/xb = 95 MeV, /i s = 17 MeV) a value Dq ~ 0.05, consistent with the measured value 0.065. 

As this discussion demonstrates, the current NA35 strange particle data are neither 
precise enough nor do they provide sufficiently complete information to allow for a similar 
analysis as done for the WA85 data above. In particular the lack of data on S~ and 2~ 
production in S-S collisions prevents us from determining the values of 7 S and /i s . In the 
absence of such information, the validity of a HG description for the NA35 S-S data can at 
the present moment not be excluded. 
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5 Implications for the fireball evolution dynamics 

So far we have extracted, with the help of thermodynamical concepts, from various sets of 
nuclear collision data the characteristic statistical properties of the source of the emitted 
particles. In the following we will comment on what our findings imply for the possible 
formation of a quark-gluon plasma, for its evolution and ultimately its dissociation into the 
final state hadrons. 

We have seen that a HG interpretation of the WA85 data on strange baryon production 
is not consistent with the large multiplicity densities observed by EMU05 which indicate 
a much larger specific entropy than implied by the HG picture. In this context we have 
discussed both a "purely thermal" interpretation of the data (which identifies the slope 
of the transverse mass spectra directly as the fireball temperature) and a "flow scenario", 
where part of the observed slope is due to a blueshift factor caused by transverse collective 
expansion of the fireball. Independent of which version of the HG picture we used, we found 
from the WA85 data a surprisingly large degree of strangeness saturation 7 S , ranging from 
75% in the purely "purely thermal" interpretation to 90% in the "flow scenario". These 
are hard to understand in a HG scenario, because the time scale for strangeness producing 
processes is usually believed to be much too slow to achieve that kind of saturation during 
the lifetime of the collision, unless a QGP is formed. 

Within the HG picture, the dynamical aspects connected with the flow seem to be re- 
quired for consistency with the freeze-out kinetics M. For this reason the "flow scenario" 
has been highly favored before. However, we have shown here that within the HG picture 
the value /i s = constrains the liberty to interpret the slope of the transverse mass spectra 
in terms of a temperature which is blueshifted by flow. The observed values /z s = and 
Hb/T = 1.63 are inconsistent with strangeness neutrality of the observed system of hadrons 
at the rather low value of T ~ 150 MeV resulting from the "flow" picture. To recover 
consistency with the strangeness neutrality of the fireball under these conditions, the HG 
equation of state would have to be very different from the conventional Hagedorn resonance 
gas which we employed here. While such a major modification of the physics of the hadron 
phase cannot at present be completely excluded, it would invalidate a large body of success- 
ful phenomenological work based on the conventional HG picture. Thus, if the experimental 
result /4 = should be confirmed via other particle ratios, in our opinion it practically 
eliminates within the framework of a conventional HG-EoS the flow interpretation of the 
spectra as a viable alternative to the purely "thermal" interpretation. 

We therefore conclude that within a HG picture the interpretation of the m^-slope of 
~ 210 MeV as the true temperature of the fireball appears to be the only possibility consistent 
with the constraint of (approximate) strangeness neutrality. This excludes the possibility of 
transverse flow, however, and thus causes consistency problems for the freeze-out process. 
In addition, it does not lead to the observed value of Dq. Therefore a HG interpretation of 
the data appears to be excluded. 

If one tries instead to interpret the observations in terms of QGP formation, one faces 
the problem to understand how the hadronization process can lead to a final state with the 
observed values for the thermodynamic parameters. We have seen that the conditions of 
strangeness neutrality and of a fixed up-down asymmetry in general lead, at the same tem- 
perature, to different values for the chemical potentials /i u , /Xd, fa i n t ne t wo phases. While 
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a vanishing value of pi a is a natural consequence of strangeness neutrality in the QGP, it is 
not at all natural in the hadron gas, where u. B — solves the strangeness neutrality condition 
only for very specific pairs of T and /zb, and thus in general it is not a value expected to arise 
in an interpretation of the experiments involving HG fireballs. On the other hand, while the 
values of the chemical potentials assumed in a hadron gas can be directly determined from 
the measured abundance ratios between the various hadronic species emitted by the gas, 
this is not automatically true for the QGP chemical potentials. How /i u , fid, and /i s inside 
the QGP are related to the abundance ratios of the finally observed hadrons depends on the 
dynamical and kinetic details of the hadronization process. 

To illustrate this problem, let us consider the limiting case of a very long-lived and slowly 
hadronizing quark-gluon plasma, where hadronization proceeds through a mixed phase pro- 
viding intimate contact between the two sub-phases. In this case chemical equilibrium not 
only within the hadronic sub-phase (with respect to processes which transform the hadrons 
among each other) and within the quark sub-phase (with respect to processes transforming 
quarks and gluons into each other), but also between the two phases (with respect to pro- 
cesses which create hadrons out of quarks and gluons or dissociate them again into quarks 
and gluons) can be maintained. As a consequence, the chemical potentials change smoothly 
from their values originally assumed in the QGP phase into the chemical equilibrium values 
of the hadron gas, and only the latter (and not the former) would be observable by studying 
particle ratios. In such a scenario, the value of the strange quark chemical potential 
extracted from a suitable combination of strange hadron abundances || would in general 
not be zero (although it started out as zero in the original QGP phase), but would be given 
by the solution to the strangeness balance constraint ( fH] ) in an equilibrated hadron gas. 

However, as demonstrated in this paper, the strange particle yields from the WA85 



collaboration [gfj indeed do point to a value of fj, s = 0, and in view of the naturalness of this 
value in the context of a QGP we are inclined to believe that this result may not be entirely 
accidental. Therefore we ask under which conditions u. s — could have survived the QGP 
hadronization process after all. 

One possibility to transfer the QGP values for the chemical potentials and temperature 
directly to the final state hadrons is to assume rapid break-up of the QGP combined with 



statistical quark recombination [|3], [13], 3T|, |32|. A necessary condition for the preservation 
of the QGP thermodynamic parameters is that sequentially produced hadrons escape suf- 
ficiently rapidly such that there is no formation of a intermediate HG phase in which the 
chemical potentials and/or momentum distributions could re-equilibrate. One can view this 
process as the peeling of the QGP surface into the vacuum — since there is no barrier, 
the surface regions of the QGP fireball are peeled off in form of rapidly escaping hadrons, 
carrying in their abundances the information about the thermal and chemical properties of 
the QGP. 

One problem with such a picture is that the system of radiated hadrons, described by 
thermal distributions with chemical potentials and temperature as given originally in the 
QGP phase, by definition is characterized by the same particle abundance ratios and by 
the same value of the specific entropy as an equilibrium HG at the same temperature. The 
radiation obtained from this quark recombination picture is in general neither strangeness 
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neutralta, nor does it carry away the correct amount of entropy per baryon (as fixed by the 
original QGP) It is well known that the entropy balance can be saved by considering the 
entropy contained in the gluons of the QGP, letting them materialize as hadrons via gluon 



fragmentation into additional quark-anti-quark pairs. In |T3| the additional quarks from 
gluon fragmentation where hadronized according to the same laws of combinatorics as the 
thermal quarks from the QGP; as a result the QGP chemical potentials were not preserved, 
but the quark chemical potentials were modified during hadronization in the direction of 
their HG equilibrium values. 

The minimal dynamical picture could consist of a primordial explosive production of the 
few high m± strange anti-baryons, followed by equilibrium hadronization of the remaining 
fireball at nearly constant entropy. To test such a picture requires the simultaneous mea- 
surement of the abundances and momentum spectra of (strange and non-strange) baryons 
and anti-baryons and of (fully identified) pions and kaons. Only with such a complete set of 
data on produced hadrons can a full understanding of the evolution of the primordial dense 
phase (QGP?) be achieved. 

Of course, our contention that the presently available data are inconsistent with a HG 
picture for the collision and rather point to a very entropy-rich source of particle emission, 
possibly a quark-gluon plasma, must appear provocative at the present stage. Our findings 
show that it is necessary to widen the investigation of strange anti-baryon ratios to higher 
and lower nuclear beam energies and to study simultaneously the particle multiplicity and 
charge flow at high and low mi, in order to confirm the now very suggestive conclusion of 
QGP formation in the collisions of S ions at 200 GeV A with heavy nuclear targets at near 
zero impact parameter. In view of the ambiguity between the HG and QGP picture, as 
far as the consistency of a vanishing strange quark chemical potential with the condition of 
strangeness neutrality is concerned, the presently available beam energies at CERN turn out 
to be somewhat unfortunate, due to the observed particular value of around 210 MeV for 
the m_i_-slope; this ambiguity could be easily resolved if /x s = were still observed in cases 
where the m^-slope (apparent temperature) is lower or higher. It is natural to expect such a 
variation in the slope parameters (combined by an inverse variation of the baryon chemical 
potential) once collisions at lower or higher beam energies are studied (for example with a 
50 GeV A Pb-beam at the CERN SPS or with 100+100 GeV A Au-beams at RHIC). We 
are also looking forward with anticipation to the forthcoming Pb-Pb collisions at about 170 
GeV A projectile energy, which will allow to check and extend the results of our analysis by 
going to the largest possible thermodynamic systems using nuclear collisions. 
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